Thermodynamical Test of Non Extensive Thermostatistics 
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An ideal mixture of parahydrogen (with nuclear spin K = 
0) and orthohydrogen (with A' = 1), in statistical weights 
1/4 and 3/4, respectively, is used as a test ground for the 
existence of non-extensivity in chemical physics. We report 
on a new bound on the non extensivity parameter q — 1 that 
characterizes generalized thermostatistics a la Tsallis. This 
bound is obtained on the basis of laboratory measurements of 
the specific heat of hydrogen. Suggestions are advanced for 
the performance of improved measurements. 
PACS: 05.30.-d, 05.70.-a, 05.70.Ce 

In the last twelve years a considerable amount of ef- 
fort has been expended investigating systems for which 
the traditional Boltzmann-Gibbs-Shannon (BGS) statis- 
tical treatment, for a variety of circumstances, seems to 
fail, al least in some respects [El. The so-called Tsallis 
non-extensive thermostatistics (TT) has played a promi- 
nent role in these endeavors TT is characterized by 
a real parameter q, such that q — I yields the standard 
BGS treatment. The first successful applications of such 
formalism was the finding of a non-divergent, physical 
distribution function for stellar polytropes j|] . A signifi- 
cant TT step was that of Boghosian. His analysis of the 
relaxation of a two dimensional pure electron plasma ^] 
yielded a g- value quite different from unity. 

Many subsequent advances, both on the basic theory 
front and in different applications, have transformed this 
non-extensive treatment into a new paradigm for statisti- 
cal mechanics TT is based upon Jaynes' MaxEnt 
approach to statistical mechanics M. Its main feature 
is that of replacing the conventional Shannon logarith- 
mic information measure S (subject to appropriate con- 
straints derived from the available a priori knowledge) by 
Tsallis' non-extensive measure (we normalize all quanti- 
ties to A: = 1) 



S, = iq- 1) 



w 



(1) 



where the set {pi}, to be obtained by recourse to Jaynes' 
variational procedure, yields the probability distribution 
associated to the relevant W microstates. Notice that for 
q — 1 one has Sq = S. 



Sq is non-extensive on account of the fact that, 



for 



a system composed of two noninteracting subsystems 

Sq{A + B)= Sq{A) + Sq{B) + (1 - q) S q{A) S q (B) , (2) 

so that q measures the degree on non-extensivity. During 
almost 10 years, the concomitant non extensive thermo- 
statistical considerations were based upon the theoretical 
framework of Curado and Tsallis 0. Nowadays, it is be- 
lieved that this framework has been superseded by a new 
one, called "the normalized one", advanced in and 
formalized in great detail by Tsallis, Mendes and Plastino 
which seems to exhibit important advantages [ p^ . 
This normalized treatment, in turn, has been consider- 
ably improved by the so-called "Optimal Lagrange Multi- 
pliers" (OLM) approach Q. The OLM treatment seems 
to provide one with a natural bridge that connects Tsal- 
lis' thermostatistics with thermodynamics [pd]-p^. The 
main difference between the Curado-Tsallis approach 0| 
and the "normalized-OLM" one is to be found in the 
manner of computing expectation values |pd]| . One has, 
for a given observable A that takes the value for the 
microstate i 



{A)q = ^p|ai, (Curado-TsaUis), 



and 



{A)q 



T,^=l Pi 
l^i=l Pi 



(normalized-OLM) . 



(3) 



(4) 



* Argentine National Research Council 



Differences in some details of the variational procedure 
distinguish the approach of Ref . |9j from the one of , 
that is to be the focus of attention in this work. 

It is quite clear that Sq provides us with an exceed- 
ingly useful information theory tool. Can we identify 
it with a thermodynamic entropy as well? For an ap- 
propriate answer an experimental determination of the 
g-value within a purely thermodynamics setting becomes 
necessary. One such a test was reported in Refs. [|l5|,0, 
where bounds to l^ — 1| were established using the cosmic 
black-body radiation. Also, present-day determinations 
of the Stefan-Boltzmann constant a put a similar con- 
straint, for which the order of magnitude is \q—l\ < 10^"*. 
Finally, bounds on the non extensivity of fermion sys- 
tems have been reported 17|, namely |<7 — 1| < 2 x 10~^. 
All these works refer to the Curado-Tsallis formalism. 
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With respect to this treatment, the OLM reformulation 
of the theory ||ll|] entails a sensible physical difference: in 
the latter, non extensivity is restricted just to Sq, while 
the energy behaves as an extensive quantity (it is non- 
extensive in the Curado-Tsallis instance). To our knowl- 
edge, no comparison of predictions arising from this OLM 
scenario with experiment has been carried out yet. 

Now, except for the measurement of Stefan-Boltzmann 
constant, all the above reported bounds on g — 1 have 
been based on observational rather than on experimen- 
tal, chemical physics data. The advantage of laboratory 
measurements in the comparison of a new theory with na- 
ture lies in the possibility of varying and controlling the 
conditions in which the data is acquired, thus ensuring 
a thorough validity test. With regards to the observa- 
tional tests, mostly of cosmological nature [p"5|-p^, it has 
been forcefully argued [|9) that they can not test the de- 
gree of non extensivity of a system, as the application of 
thermodynamic arguments is of a strictly local character. 
This entails that limits on the non-extensivity, in a large, 
cosmological scale, can not be thereby obtained. 

In order to test non extensivity, a (laboratory) physical 
system composed of two different weakly interacting sub- 
systems should be used. It is well known that H2 gas is 
a mixture of two different substances : parahydro- 

gen (with nuclear spin K = 0) and orthohydrogen (with 
K = 1), in statistical weights 1/4 and 3/4, respectively. 
Since there exist accurate (about 1%) measurements of 
both its specific heat and relevant spectroscopic data, this 
ideal mixture is a very good test ground for the existence 
non-extensivity in chemical physics. In this communica- 
tion a new bound on the non extensivity parameter q — 1, 
based on laboratory measurements of the specific heat of 
hydrogen, is reported, and suggestions are advanced for 
the performance of improved measurements. 

Since the proposal of Giauque , it has been custom- 
ary to compute the thermodynamic properties of simple 
substances, such as hydrogen, from spectroscopic data. 
Indeed, when accurate spectroscopic data are available, 
this method gives much greater precision than direct 
thermal measurements. On the other hand, the method 
rests on the validity of BGS statistics as the correct de- 
scription of nature. Giauque himself computed the ther- 
modynamic properties of the II2 molecule [ p3| and found 
very good agreement with the measured data [p^-p6|. 
The Giauque proposal has been very successful, to the 
point that no new direct measurement of the specific heat 
of II2 molecule has been reported since 1930! 

At this point we start our non-extensive considerations. 
Since the experimental thermodynamic value of q can 
not be expected to appreciably differ from unity (other- 
wise, non-extensivity would have been discovered decades 
ago!), a perturbative treatment around q = 1 will be un- 
dertaken. In order to compute the specific heat of a sub- 
stance in the standard Curado-Tsallis (CT) formalism, 
the point of departure is given by the Tsallis partition 
function 



where 



+ 



(5) 



1/T with T the 



temperature and {ej} is the set of energy levels. We 
compute then the Tsallis free energy 

1 



q 



(6) 



(x^ * — 1)/(1 — g), and finally, apply the 
usual thermodynamic relations 



where In^ x 



c(CT) ^ '^^1 
« dT ' 



« dT ' 



(7) 



An straightforward expansion of equations (|^) and (|^) 
around q—1 yields 



Ff^^=F-\{q-l) 



(8) 



where U and F are the BGS {q = \) internal and free 
energies respectively and C is the BGS specific heat. 

Additionally, for an out of equilibrium (/, 1 — /) mix- 
ture of two substances (A,B), the specific heat satisfies 
the analog of equation 



Cf ^)(A + B)^ Cf^\A) + Cf^\B) + (1 - q)C 
Cab 



AB. (9) 

^) (A)5f ^) [B) + ^) (B)5f ^) [A)] . (10) 



The calculation of the specific heat according to the 
OLM treatment is of a somewhat more subtle nature 
because the thermodynamic relations does not keep the 
usual form in terms of Tsallis' entropy One starts 
from the definition 



dT ~ ^ d/3' 
where the internal energy Uq can be written as |ll[ 



Uq^ 



(11) 



(12) 



In this equation, the fj — 1 — (1 — q)J3^j arc the 
so-called configurational characteristics |ll[ where Aj — 
tj — Uq. Taking the pertinent derivatives, one arrives to 



Ca 



9/32 l-qf}{^,f-^) 



(13) 



which in the q ^ 1 limit goes into the BGS specific heat. 
No corrections for the fact that we deal with a mixture 
are necessary, since the internal energy is additive in the 
OLM approach fll 



U, 



OLM^ 



{A + B) = t/"^^(A) + f/"^^(B) 



(14) 
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The specific heat of parahydrogen (orthohydrogen) can 
now be written in the fashion 



with 



C, = C + (q - l)AC\ 



dC 



(15) 



AC^'^^' =2piFC + Ul3 ,^ 

V dp 

dC 



(16) 



for the Curado-Tsalhs version of the formahsm, and, with 



A| = (e, - [/), 



AC" 



OLM 



( \ (Af 



(17) 



for its OLM counterpart. Notice that all quantities en- 
tering AC are evaluated within the framework of the 
standard, BGS statistics. For the case of H2 one has now 



AC(^^)(H2) = ^ACf ^' + f ACf ^) 



AC 



OLM 



(H2 



A /^OLM 



(18) 
(19) 



where Co,Cp (6*0, 5'p) are the BGS specific heats (en- 
tropies) of ortho- and parahydro gen , respectively. 

The structure of AC in Eq. (Hq) exhibits in a quite 
clear fashion i) a characteristic non-extensive behavior 
and ii) a dependence on the reference point for the energy. 
In Eq. (^^ , on the other hand, the non-extensivity is not 
apparent! It is hidden in the higher order momenta of the 
energy. Also, the troublesome dependence on the origin 
of the energy scale just vanishes. 

Table || displays the accurate data used in this work. 
The theoretical BGS thermodynamic functions were 
computed from spectroscopic data using the usual ex- 
pressions. The constants for the rotational levels were 
taken from Ref. |^^, since they are accurate enough to 
get a good representation of experimental data. The O- 
C differences between BGS theory and experiment (Fig. 
|l|) were used to fit the Tsallis parameter (g — 1) with a 
simple least squares procedure. The ensuing results are 
shown in the first line of table ||. 

These results, however, may be contaminated by sys- 
tematic errors. The largest deviations from the BGS spe- 
cific heat come from the data in reference ||2^. Assum- 
ing that in this data there is a small deviation from the 
3/4 — 1/4 mixture, the O - C differences would adopt the 
appearance 



ACofc3 = (<Z-l)AC + e(Cp-Co) 



(20) 



where e is an auxiliary parameter. In such a case one 
obtains the results of lines 2 and 3 of table ||, which are 
consistent with zero at the 95% CL. 

This suggest the existence of systematic errors in the 
data. 

As a final result, we can state that the Tsallis param- 
eter is consistent with zero at the 95% confidence limit: 



\q-\\cT < 5 X 10"* (95%C.L.) 
1 loLM < 1.2 X 10-3 (95%C.L.). 



(21) 
(22) 



The above equations constitute the main result of this 
communication . 

Our results show that the BGS statistics is, among 
the set of theories parameterized by q — 1, the one that 
describes the available data with an accuracy better than 
10"'^. Within this level of accuracy, there is no evidence 
on a non-extensive behavior in this simple system. 

With better experiments, the accuracy of our results 
( ^ ) and ( p^ ) could be improved, in principle, by several 
orders of magnitude. Indeed, spectroscopic data have a 
potential accuracy of 10~^, while an accuracy of 10""* 
should be easily attained with the use of modern exper- 
imental techniques. In such a way a stricter test for the 
validity of the BGS statistics (or, alternatively, for the 
existence of non-extensivity) would become available for 
this system. 
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FIG. 1. 0-C residuals of the specific 
sped to the Bolztmann-Gibbs-Shannon 



heat of H2, 
theory. 



r(K) 



c 



Ref T (K) 



C 



Rcf T (K) 



C 



Ref 



35.00 1.505±0.005 |24| 94.12 1.688±0.012 
40.00 1.505±0.010 1241 100.00 1.727±0.046 
45.00 1.515±0.008 1241 100.42 1.746±0.012 
50.00 1.520±0.018 1241 106.64 1.777±0.012 
60.00 1.510±0.018 1241 110.00 1.828±0.030 
65.00 1.535±0.004 i24| 112.52 1.819±0.012 
70.00 1.565±0.007 1241 118.63 1.847±0.004 
80.00 1.586±0.004 1241 135.71 1.967±0.004 
81.12 1.601±0.004 \2e\ 136.62 1.957±0.013 
82.00 1.611±0.007 1241 142.89 1.975±0.014 
85.00 1.621±0.012 1241 145.64 2.041±0.005 
90.00 1.646±0.010 || 148.98 2.041±0.014 

TABLE I. Experimental data on the H2 specific heat. The 
columns display, respectively, the temperature in K, the molar 
specific heat in units of R with our adopted standard devia- 
tion, and the corresponding bibliographic reference. 
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Parameter CT Version OLM Version 



q-i 


-0.21 ±0.05 


0.89 ±0.23 




-0.11 ±0.07 


0.52 ±0.57 


e 


1.0 ±3.7 


2.9 ±4.1 



TABLE IL Results of the analysis performed in this work. 
The columns of the table display the adjusted parameters' 
names and their values (in units of 10~^) for the two versions, 
CT and OLM, of the theory. The first line lists a raw fit to the 
residuals, while the second and third lines display a fit that 
assumes a small departure from the 1/4 — 3/4 proportions of 
the mixture. 
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